Correlated Fluctuations between Luminescence and Ionization in Liquid Xenon 
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The ionization of liquefied noble gases by radiation is known to be accompanied by fluctuations much 
larger than predicted by Poisson statistics. We have studied the fluctuations of both scintillation and 
ionization in liquid xenon and have measured, for the first time, a strong anti-correlation between 
the two at a microscopic level, with coefficient —0.80 < p ep < —0.60. This provides direct experi- 
mental evidence that electron-ion recombination is partially responsible for the anomalously large 
fluctuations and at the same time allows substantial improvement of calorimetric energy resolution. 
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The measurement of ionizing radiation in a liquiefied 
noble gas 1 such as Ar, Kr, or Xe can be characterized 
by two parameters: W e , the mean energy required to 
create a free electron-ion pair and the Fano factor, F ei 
which parameterizes the fluctuations in the number of 
ion pairs. 2 The Fano factor is defined by 



a e = y/F e N e 



(1) 



where cr 2 is the variance of the charge expressed in units 
of the electron charge, e, and N e is the number of ion 
pairs. F e = 1 corresponds to Poisson statistics. Fano 
originally predicted that the charge fluctuations would 
be sub-Poissonian, F e < 1, because the individual ion- 
pair creation processes are not independent once the ad- 
ditional constraint E — N e W e is included, where E is 
the energy deposited by the incident particle. The Fano 
factor for gaseous noble elements is relatively well under- 
stood, where for example the theory for argon predicts 
F e — 0.16 while experiment yields F e = 0.20. 3 

The Fano factor for the liquid phase has been far more 
difficult to understand, even though the values of W e are 
similar to those of the gas. The pioneering theoretical 
work of Doke 4 in 1976 predicted F e w 0.05 for liquid Xe 
(LXe). Experiments, however, showed F e > 20, that 
is, a variance 20 times worse than Poissonian and 400 
times worse than predicted by Fano's original argument. 
This discrepancy has not only raised important issues 
for the understanding of the phenomena of energy loss 
and conversion, but it is also of interest for experimental 
physics since this large Fano factor limits the resolution 
of calorimeters used in nuclear and particle physics. 

Luminescence light provides a second process, comple- 
mentary to ionization, with which to study energy loss 
and conversion phenomena. The properties of lumines- 
cence light emitted by ionizing radiation (scintillation) 



in liquid Ar, Kr, and Xe detectors have been extensively 
studied and have been exploited in calorimetry 1 ' 5 . As 
described elsewhere 6 , scintillation photons are produced 
by the relaxation of a xenon excimer, Xej — > 2Xe + 7. 
The excimer is formed in two ways; 1) direct excitation 
by an energetic particle, 2) recombination of an electron 
with a xenon molecular ion, Xe^ + e — ► Xei, (recombina- 
tion luminescence) 7 . The electron-ion recombination rate 
depends upon the ionization density and the applied elec- 
tric field. The quantities W p and F p can be defined in 
analogy with the ionization case, where W p is the mean 
energy absorbed per emitted photon and F p parameter- 
izes their variance. Fano's theory was not intended to 
treat photon fluctuations and in this case F p is treated 
as a phcnomcnological parameter. 

We have simultaneously studied both ionization and 
scintillation in LXe using a special geometry ionization 
chamber equipped with a VUV photomultiplier tube 
(PMT) capable of detecting the ~ 175 nm scintilla- 
tion light from the xenon. The chamber, shown in Fig- 
ure 1, has a pancake geometry with a total volume of 
~ 1.5 liters. It is built out of ultra- high- vacuum com- 
patible parts, so that it can be baked and treated with 
hot xenon before filling to obtain good LXe purity. The 
chamber is cooled by immersion in a bath of HFE7100 8 
that is in turn cooled by liquid nitrogen. The xenon is 
purified by an Oxysorb 9 cartridge during transfer from 
the storage cylinder to the chamber. Ionization is col- 
lected in a fiducial volume with radius ~ 10 mm and 
height ~ 6 mm defined by a small grid and a large cath- 
ode plane. The grid and cathode are made of an elec- 
troformed nickel mesh with 90% optical transparency. 
Drifted electrons are collected on a low-capacitance an- 
ode surrounded by a field defining ring held at the same 
ground potential. An electron transparency of 100% is 
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FIG. 1. Schematic drawing of the liquid xenon detector. 
As described in the text, the chamber is cooled by immersion 
in a refrigerant bath. The bismuth source is in reality much 
smaller than indicated in the figure. 



obtained by maintaining the anode-to-grid field, (Eag), 
double in value with respect to the drift field, (Ekg) 10 , 
for all the measurements reported here. All data was col- 
lected at a LXe temperature of 1 69. 2 ±0.5 K and pressure 
of 990^5q Torr. The temperature gradients in the cham- 
ber were measured to be less than ±0.2 K. 

The electron signal is detected by a low noise charge 
sensitive preamplifier 11 with a cold input FET placed on 
the feedthrough of the chamber. After subsequent am- 
plification and filtering, the charge signal f(t) is recorded 
on a transient digitizer and fit to a model assuming drift 
to the anode of a gaussian ionization distribution repre- 
sented by the function 



f(t) = A 



I + " / ( ^ 



+ a 



(2) 



where At = t-t and the term a = De~~ + 1 00 + 0.019* 
accounts for pile-up, ADC offset and a small integral 
non-linearity. We simultaneously fit for the amplitude A, 
the effective drift-time t$, the rise-time parameter C, the 
pile-up offset D and the integrator fall-time constant r. 
The fit fall-time is ~ 50 fjs for all the data, as expected. 
The parameter C includes information about the multi- 
plicity of the ionization event and its values range from 
0.78 /US for a field of 4 kV/cm to 1.04 fis for 0.2 kV/cm. 
The decrease of this parameter with the increase of the 
electric field is consistent with the expected change in 
drift velocity. 

The readout is calibrated in units of electron charge 
by injecting a signal into the input FET through a cali- 
brated capacitor. The noise of the preamplifier was mea- 
sured with a test signal and found to follow a gaussian 
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FIG. 2. The two-dimensional scintillation and ionization 
spectra recorded at drift field Ekg = 4kV/cm. The two "is- 
lands" with negative correlation coefficient correspond to the 
two 7 lines from the 207 Bi source and their satellite internal 
conversion peaks. The axes are calibrated in terms of abso- 
lute numbers of elementary excitations (ionization electrons 
and photo-clcctrons in the PMT). 



distribution with width <tn = 381e. 

The scintillation light is detected through the cathode 
grid by the 3-inch VUV PMT 12 with a solid angle that 
varies from 21.5% to 26% of 4-7T for signals produced in 
the region between cathode and grid active for ionization 
collection. The fraction of photoelectrons extracted from 
the PMT photocathode per scintillation photon produced 
in the LXe varies between 3.0% and 3.6% considering 
mesh and quartz window transparencies and photocath- 
ode quantum efficiency (22% including transparency of 
the PMT window). The PMT signal is read out using a 
12 bit charge ADC calibrated by measuring the single- 
photon peak, which is well resolved from noise. 

Data acquisition (DAQ) is triggered when there is a 
delayed coincidence between PMT and charge signals. 
The delay is required to account for the electron drift- 
time and varies with electric field. Once a coincidence 
is established, DAQ is inhibited until the event has been 
processed and the instruments cleared. About 25% of 
the collected data is removed off-line by analysis cuts 
that reject events that don't fit model (2). These "bad 
events" are mostly due to multiple ionization deposits 
from Compton scattering, deposits at the edge of the 
ionization fiducial region, and cosmic ray background. 

The LXe is excited by the 570 and f064 keV 7's 
(and their associated internal conversion electrons and 
X-rays) emitted by a 207 Bi source. The source is pre- 
pared as a monoatomic layer electro-plated on a ~ 3/im 
diameter carbon fiber that is then threaded through the 
cathode grid. The thin source minimizes the effects 
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of self-shadowing of the scintillation light and is essen- 
tial to obtain acceptable scintillation resolution. The 
data presented here were collected at drift fields Ekg of 
0.2kV/cm, 0.5kV/cm, lkV/cm, 2kV/cm, and 4kV/cm. 

A two-dimensional scatter plot of ionization versus 
scintillation signals obtained at 4 kV/cm is shown in Fig- 
ure 2. The two peaks in the ionization channel can be 
fit to a gaussian function plus a first order polynomial 
function. The energy resolution obtained from the ion- 
ization data, once the electronics noise is subtracted in 
quadrature, is in good quantitative agreement with the 
results of other authors, as shown for the 570 keV peak 
in Figure 3. 
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FIG. 3. Ionization resolution for 570 keV 7s in LXe when fit 
to a one-dimensional gaussian. The black squares show the 
resolution observed by this experiment for different electric 
fields. The other symbols represent the resolutions observed 
with LXe ionization chambers by other authors , in good 
agreement with our data. 

As shown in Figure 2, the scintillation spectrum has 
substantially poorer resolution and a larger continuum 
than the ionization spectrum. Differences between the 
two channels arise from the fact that the VUV photon 
acceptance depends somewhat upon the position of the 
energy deposited, as mentioned above, and from the dif- 
ferent fiducial volumes. The resolution in the scintilla- 
tion channel, as extracted by a fit with a gaussian plus a 
first degree polynomial, is a — (16.5 ± 2.5)% at 570 keV, 
where the error refers to the range of values observed at 
different fields. 

The observed ionization and scintillation spectra are 
qualitatively well reproduced in both peak positions and 
strengths by a detailed GEANT 4 14 simulation which in- 
cludes the processes of internal capture, X-ray and Auger 
electron emission. An empirical gaussian smearing is 
used to match the GEANT 4 output with the data. How- 



ever, in the absence of a microscopic model for the pro- 
duction of scintillation photons in LXe, the simulation 
could not be used to quantitatively fit the two dimen- 
sional spectrum. 
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FIG. 4. Macroscopic anti-correlation between scintillation 
light and ionization as a function of the electric field. Similar 
results have been obtained before by other workers. 

Previous works 7 have reported evidence for what we 
call here "macroscopic" anti-correlation between scintil- 
lation light and ionization as a function of the electric 
field applied to the LXe. This effect, illustrated for our 
data in Figure 4, is thought to be due to the fact that 
scintillation occurs as ionization rccombincs. Different 
values of the electric field modify such recombination, 
changing the proportion between light and ionization but 
leaving their sum constant. The constancy of the sum is 
well reproduced by data, as shown, in our case, by Fig- 
ure 4. We note that the absolute value of the sum does 
not directly correspond to the total energy deposited, as 
part of the energy is lost in other channels. We refer to 
this anti-correlation as "macroscopic" because is it an av- 
erage behavior of the LXe and it is obtained by changing 
an external, macroscopic, parameter (the electric field). 

In Figure 2, clear "microscopic" anti-correlation is ev- 
ident in each of the two "islands" corresponding to the 
570 keV and 1064 keV 7 lines. This effect, simply referred 
to as anti-correlation below, indicates that energy sharing 
between scintillation and ionization fluctuates event by 
event even though the sum of the two remains constant, 
in analogy with the macroscopic correlations described 
above. The microscopic nature of the anti-correlation 
may help in understanding the large Fano factor men- 
tioned above. It also shows the way to dramatic improve- 
ments in the energy resolution of particle detectors, since 
the fluctuations clearly affect scintillation and ionization 
in opposite ways. Anti-correlations similar to the ones 
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shown in Figure 2 are obtained at all fields and were not 
previously observed in LXe, although they were shown, 
without a detailed discussion, for heavy ions in liquid 
Argon, by Ref. 15 . 

The anti-correlation can be formally extracted from 
the data by fitting each of the peaks to a 2-dimcnsional 
gaussian distribution 16 



2p ep AN e AN 



G(N e , N p ) = He 



(3) 



1 



p\p is a normalization con- 



where H = l/2ircr' e (7p^ 

stant, p ep — V e p/a' e a p is the correlation coefficient be- 
tween ionization and scintillation signals, V ep = ((N e — 
N e )(N p - N p ))js the covariance, AA e = N e - N e and 
AN P = N p — N p . While the instrumental noise on the 
scintillation channel is negligible and hence ignored, the 
variance for the ionization channel, a 12 = a 2 + a 2 ^, in- 
cludes the term, un — 381 e, that accounts for the read- 
out noise mentioned above. This is a relatively small cor- 
rection (e.g. 12% at 4 kV/cm for the 570 keV peak). Un- 
like the one-dimensional fits used for ionization and scin- 
tillation data separately, the 2-dimensional fitting func- 
tion does not include a background term. A two-step pro- 
cedure is used to first find an elliptical region of full width 
2.35cr around each of the two 7-ray peaks and then per- 
form the final fit. The limited range of the fit suppresses 
the K-shell internal capture electrons and Compton scat- 
tering events. We found that the two-dimensional gaus- 
sian fits tend to produce larger variances with respect 
to the one-dimensional fits, providing a conservative es- 
timate of the real resolutions. 

The correlation coefficients found by the fit are plot- 
ted in Figure 5 as a function of drift field. A clear anti- 
correlation exists at all fields for both the 570 keV and 
the 1064 keV lines. The errors shown in the Figure are 
systematic and are estimated as the full range of the vari- 
ation of p ep when the different analysis cuts and fit in- 
tervals are changed. Statistical errors are negligible. The 
stability of the correlation with respect to different exper- 
imental circumstances (such as position in the chamber, 
size of the chamber and type of radiation) can be tested 
by repeating the fits for different subset of the data, each 
presenting a different range of effective drift-times, as 
derived from Eq. 2. The effective drift-time distribution, 
shown for 4 kV/cm in the inset of Figure 6, allows for 
the selection of events originating at different spatial po- 
sitions along the drift field. The peak at large drift-times 
primarily corresponds to electrons that lose their energy 
and stop in the vicinity of the cathode plane where the 
source is. As shown in Figure 6 the correlation coefficient 
does not depend significantly on different effective drift 
times. 

By diagonalizing the covariance matrix via a unitary 
rotation through angle 9 — 2p e pa e a p / (a 2 — a 2 ), the cor- 
related two-dimensional gaussian can be reduced to two 
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FIG. 5. Correlation coefficients between ionization and 
scintillation signals obtained from the 2-dimensional fit dis- 
cussed in the text. The two sets of points refer to the coeffi- 
cients found for the 570 and the 1064 keV peaks. The errors 
shown are from systematics in the analysis, while statistical 
errors are negligible. 
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uncorrelated gaussians with variances a 2 
This rotation is mathematically equivalent to measur- 
ing the variances along the major and minor axes of the 
two-dimensional ellipses. Figure 7 shows the ionization 
resolutions from the 2-dimensional fit and the minimum 
resolutions in the frame rotated by 8 for the 570 keV 
gamma ray peak at all drift fields. 

The resolution improvement observed when correlat- 
ing the ionization and scintillation signals is significantly 
better than the uncorrelated addition of the two spec- 
tra, where cr p+e = {o~ 2 + a^ 2 )^ 1 / 2 . At 4 kV/cm 
3.7%. When using the correlated nature of the 



3.0%, a factor of 



Up+e 

two channels the resolution ay, 
1.24 improvement. The improvement is more dramatic 
at the 0.2 kV/cm field where the ratio is 9.7/5.4 = 1.8. 
While we found that the absolute values of the resolution 
change somewhat for different fitting methods and at dif- 
ferent stages of the iterative procedure, the improvement 
factor is very stable. We note here that the main results 
of this work are to show that a microscopic correlation 
exists and it can be consistently used to improve the reso- 
lution. The ultimate resolution will have to be measured 
at a later stage with a detector of larger volume, better 
VUV light collection and, possibly, a different source and 
it is expected to be even better than the resolution <j m in 
reported in Figure 7. 

While the nature of the large Fano factor in lique- 
fied noble gases is still not well understood, the different 
models considered generally involve fluctuations in the 
density of the energy deposited 17 . Such fluctuations can 
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FIG. 6. Correlation coefficient at all fields for different sub- 
set of events at a 4kV/cm field. The complete data set for the 
570 keV region, is divided in four regions, according to their 
drift-times, as illustrated by the effective drift-time distribu- 
tion in the inset. The peak to the right of this distribution 
is mostly due to electrons that stop near the source on the 
cathode. For each field the cuts defining different regions are 
set in such a way as to select equivalent populations of events. 
The magnitude of the correlation coefficients and their behav- 
ior as function of the electric field do not depend significantly 
upon the energy deposition location. 



arise from density variations in the medium or from the 
small number of delta electrons with energies between 1 
and 30 keV that are produced by Rutherford scattering 
of atomic electrons in the liquid. When these electrons 
come to a stop they produce a dense cluster of electron- 
ion pairs whose recombination rate is relatively high. In 
both cases more recombination occurs where the ioniza- 
tion density is higher. Hence the statistics relevant for 
the fluctuations in the liquid is that of the (small) num- 
ber of high ionization density regions (either the number 
of delta-electrons or the number of high density regions). 
In a sense, this would mean that the Fano factors com- 
puted in Eq. 1 are too small because the "wrong" N is 
used. 

As first suggested in Ref. 18 the recombination hypoth- 
esis above implies that the fluctuations of ionization 
and scintillation should be anti-correlated, since the high 
ionization-density regions which produce less free charge 
should emit a correspondingly larger amount of scintilla- 
tion light. A linear combination of ionization and scin- 
tillation should therefore have smaller fluctuations than 
either alone. Note that since the Fano factors F e and 
F p are typically far in excess of their Poissonian val- 
ues, the improvement should in principle be substantial 
and is not limited by the factor of \[2 one would ex- 
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FIG. 7. Resolutions for the ionization channel and for the 
small correlated combinations using a fit to a two-dimensional 
correlated gaussian. The data shown is for the 570 keV 
gamma ray peak. 



pect from the combination of two uncorrelated signals of 
equal resolution. The data presented here provides direct 
experimental evidence of this fluctuation and the associ- 
ated recombination by observing a clear microscopic anti- 
correlation between ionization and scintillation. The cor- 
relation, p ep , and the magnitude of the variances, <T^ ax 
and cr^ lin , provide input to improve the model in terms of 
observables that are more directly related to the micro- 
scopic physics. It is interesting to note that from Figure 5 
it appears that the magnitude of the correlation p ep de- 
creases for increasing electric field (at least at low fields) . 
This behavior is possibly connected to the dependence of 
the intensity of the recombination luminescence with the 
electric field, as expected in all the models discussed. 

In conclusion, we have observed a clear microscopic 
anti-correlation between ionization and scintillation in 
liquid xenon. We have also measured the correlation co- 
efficient between these two channels to be —0.8 < p ep < 
—0.6. This observation has two important consequences. 
First, it provides direct experimental evidence that re- 
combination is a substantial source of the anomalously 
large Fano factors in the high density liquid state. Sec- 
ond, it provides a method for improving the resolution of 
liquefied noble gas calorimeters. 
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on the subject, M. Moe for a careful reading of the 
manuscript and to R.G.H. Robertson for having stimu- 
lated the interest on this topic. This work was supported 
by US-DoE grant DE-FG03-90ER40569-A019 and by a 
Terman Fellowship of the Stanford University. 
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